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Quantum light propagation in longitudinally inhomogeneous media 
as a spatial Lewis-Ermakov physical invariance 
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We study the propagation of quantum states of light in separable longitudinally inhomogeneous 
media. By means of the usual quantization approach this kind of media would lead to the unphys¬ 
ical result of quantum noise squeezing. This problem is solved by means of generalized canonical 
transformations in a comoving frame. Under these transformations the generator of propagation is 
a physical Lewis-Ermakov invariant in space which is quantized and, accordingly, a propagator con¬ 
sistent with experiments is obtained. Finally, we show that the net effect produced by propagation 
in these media is a quantum Gouy’s phase with application in quantum interferometry. 
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I. INTRODUCTION 

In the last years great attention has been drawn to the 
time dependent quantum harmonic oscillator. This prob¬ 
lem has been thoroughly studied both from mechanical 
and electromagnetic points of view, showing 

squeezing effects. The approach most used in the tackle 
of this sort of problems is the use of quantum mechani¬ 
cal invariant operators, in particular the Lewis-Ermakov 
invariant Q. 

In the field of guided propagation of light, the anal¬ 
ogous problem is the study of longitudinally inhomoge¬ 
neous {LI) media. An invariant approach was exten¬ 
sively used in Q and more recently in Q in the study 
of graded index media but, in both cases, they studied 
classical problems by means of a quantum-theoretical ap¬ 
proach. Quantum propagation problems have been dealt 
with in homogeneous media, as shown in Q and refer¬ 
ences therein. These studies provided the background of 
the quantum theory of light propagation showing that the 
operator which describes correctly the quantum spatial 
propagation along an arbitrary direction z is the Momen¬ 
tum operator M., since the Hamiltonian approach fails 
for problems like dispersive media, counterpropagation 
and LI media as well Q. A first app roximation to LI 
media were carried out by Abram [l0| and Glauber and 
Lewenstein [llj| . where single optical discontinuities be¬ 
tween homogeneous media were analyzed. In these stud¬ 
ies they showed the different physical behaviour the fields 
experience from the analogous mathematical problem of 
the time dependent quantum harmonic oscillator, prov¬ 
ing that the field quadrature noise does not exhibit real 
squeezing, in agreement with experiments. 

As far as we know, there is not any work dealing with 
quantum states of light propagating in LI media in such a 
way that field quadratures do not exhibit real squeezing. 
So, our purpose is to give a phenomenological approach 


to this problem in waveguides. To this end, we will start 
obtaining the classical solutions which help us to find a 
reference frame which continuously performs a physical 
variable change in LI media and eliminates the virtual 
squeezing, leading us to a proper propagation genera¬ 
tor, a quantum Lewis-Ermakov type Momentum opera¬ 
tor Ad. Erom this operator we will obtain Eock states 
in the optical-field strength (OES) representation £ and 
derive the propagator, proving the lack of squeezing in 
this representation and the arising of a quantum Gouy’s 
phase. Eurthermore, we will present the particular case 
of propagation of a gaussian quantum state in a cosine- 
type LI medium, where will be shown that the net effect 
of this kind of media on quantum states is the genera¬ 
tion of a quantum Guoy’s phase dependent on features 
of both the media and the input quantum state. 


II. CLASSICAL ANALYSIS OF PROPAGATION 
IN LONGITUDINALLY INHOMOGENEOUS 
MEDIA 


Our aim is to study the propagation of waveguided 
modes of quantum light in dispersion-free and non¬ 
magnetic media with separable inhomogeneous refractive 
index in an arbitrary direction of propagation z, given by 


n'^{x, y, z) = nl f{x, y) + An^ ^^(z), (1) 

where the longitudinal h{z) and transversal f{x,y) parts 
of the index are completely independent and uq and An 
are constants. We focus on the separable index prob¬ 
lem as it does not show coupling and therefore radiation 
modes (losses). 

Erom Maxwell equations, it is easy to show that the 
electric field E{x,y, z,t) obeys the vectorial wave equa¬ 
tion [l^ : 


d'^Et 
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with E = {Ex, Ey, Ez) = {Et,Ez). Let us consider 
monochromatic guided ID vector modes with frequency 
Wo- represented by vector held solutions with the follow¬ 
ing factorable complex amplitudes: 

Et{x, y, z,t) = Y^ q^c{z)itAx, v) (3) 

(7 

where we have used a for simplicity standing for 
the modal numbers v, y in each transverse direc¬ 
tion, z-dependent complex coefficients Qaciz) fulhlling 
X^o-l9o-c('Z)P = 1; and electric normalized transverse 
complex amplitudes $,f^{x,y) corresponding to quasi- 
TE (or quasi-TM) modes fulhlling a quasi-complete or¬ 
thonormalization condition Q, which belong to the ho¬ 
mogeneous part of the refractive index and satisfy: 

'^t^ta+f^onlf{x,y)ita+ 

(4) 

Vt{i,,Vt{\Tinlf{x,y))=Pli,,, 

with Pt the transverse propagation constant. Solutions 
of this equation give us the invariant transverse modal 
structure of the held. Applying equations ([T]), ([3]) and 
O into we obtain: 


with Poa = /3cr(0). Equation ([8]) is an Ermakov-Pinney 
equation with solutions given by [Tsj l: 


Po{z) = [{Poo Ua{z)f -b vl{z)Y/'^, (10) 

p,(0) = l, pU0) = 0, (11) 

and where u„ and are linearly independent functions 
that satisfy equation (O and have the following initial 
conditions and Wronskian: 

n^(0)=<(0) = 0, (12) 

uU0)=?;.(0) = l, (13) 

Wo =u'^Va - v'^Ua = 1. (14) 


So, the classical wave amplitude changes continously dur¬ 
ing propagation by a factor whereas the modal prop¬ 
agation constant is given by Pa = 0'^ via equation (|9]) 
where 9a represents the total phase accumulated in the 
propagation. 

III. QUANTIZATION IN LONGITUDINALLY 
INHOMOGENEOUS MEDIA 


+/3^(^:)9<t = 0, (5) 

where qa = {qac + <lac)l‘^ stands for the real electric held 
coefficients, Pa{z) = Pta + kQ^rPh?{z) is the local prop¬ 
agation constant of the cr-mode and where we have used 
the approximation E^ <C , Ey. It is important to out¬ 
line that Ez = 0 in the case of TE modes, that is, equa¬ 
tion ([5|) is exact for such modes. This propagation equa¬ 
tion clearly suggests a local spatial harmonic oscillator 
and therefore it can be directly derived from spatial-type 
Hamilton equations where the Hamilonian is substituted 
by the Momentum, since it is the generator of spatial 
translations Em, given by: 

Ma = ]^[pI +Pl{z)ql], (6) 

with Pa = q'a and prime stands for z-derivative. This 
result is analogous to that obtained in Q where time- 
dependent linear media was studied. The classical Mo¬ 
mentum (O is equivalent to the Hamiltonian of a time- 
dependent harmonic oscillator, with P{z) playing the role 
of w(t) 0. Likewise, the solution of equation ([5|) is easily 
obtained via the use of the complex electric field in 
the following way: 

qac{z) = Pa e*®‘'gcrc(0), (7) 


where pa and 9a are real functions obtained by solving: 



dz'^ 


I a2f \ Poa 

+ Pa{z)Po = -Z-, 

pI 

d9a _ Poa 

dz pI ’ 


( 8 ) 

(9) 


The next step is to quantize the classical Momentum 
([S]). But before carrying out this step, some considera¬ 
tions have to be taken into account. Direct quantization 
of the classsical fields qa and Pa would lead to a quan¬ 
tized z-dependent Momentum analogous to the Hamilto¬ 
nian for a time-dependent harmonic oscillator and, there¬ 
fore, following the usual steps, we would have propa^- 
tion equations leading to quadrature noise squeezing [ij- 
1^ . But in the sudy of propagation in LI dielectric media, 
this approach does not provide consistent results. As was 
pointed out by Abram in his seminal paper about quan¬ 
tization of light in dielectric media |10| , when quantum 
states propagating in a dielectric are represented in the 
basis of free-space photons, they seem to be squeezed, but 
inside a dielectric there is no experiment that can detect 
free-space photons. This happens because if we wish to 
detect photons inside the medium the fields would ex¬ 
periment an effective refractive index equivalent to the 
squeezing parameter and, therefore, a scale change is nec¬ 
essary. In the same spirit an interesting discussion about 
quantization in a dielectric was carried out by Glauber 
and Lewestein in El- In this study is stressed that the 
measurement of quantum states of light is based on pho¬ 
toabsorption processes carried out in the basis which di¬ 
agonalizes the Hamiltonian in every medium. This is so 
because the photocount distribution is given by a nor¬ 
mal ordered correlation product in the local basis, result 
of the physical property of the photoabsorption process 
that the energy of the field decreases when a photon is 
absorbed. But in any other basis this is not true, because 
of the mixing of absorption and emission operators, and 
therefore of frequencies, lacking the normal ordering in 
the correlation product. Thereby the state in these bases 
can not be measured and the photons are so-called virtual 
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and, accordingly, the squeezing is virtual as well. This 
statement is closely related with Abram’s idea and has 
to be applied to LI media as it is the continuous limit of 
single discontinuities. 

Hence, we look for some transformation which takes 
the classical Momentum to be a constant of motion 
suitable to be quantized and simultaneously get rid of 
the virtual squeezing in a continuous way. Since carry¬ 
ing out the usual approach CHI squeezing proportional 
to the function p{z) would appear, following the idea of 
a scale change in the discontinuity C3,[iil we propose 
the following generalized canonical transformations from 
a generating function G 2 {q,P) = Pq/p and z transfor¬ 
mation 0 : 

pZ 

Qa = —, P<7 = PaPa, Sa = { pl{z) dz)~^; (15) 

Per Jo 

and the next gauge transformations related to the new 

s-frame: 

Qa=Qa, Pa=Pa-—Qa, (16) 

Pa 

with the dot standing for an s-derivative. Applying these 
relations into we directly obtain the following Mo¬ 
mentum: 

MAQ,P,s) = l[P^is) + PlQlis)], (17) 

where we have used the auxiliar equation (l 8 |). It is re¬ 
markable that these new phase space and longitudinal co¬ 
ordinates (Q, P, s) are analogous to the comoving coordi¬ 
nates and the conformal (or arc-parameter) time used in 
cosmology and, likewise, the canonical momentum trans¬ 
formation (ITSl) is related to the Hubble’s law CB- After 
these transformations, the continuous change in the value 
of the harmonic potential ap pea rs as constant in the new 
comoving frame of reference jl9j. So, in this frame, there 
are not scale changes differentially and the medium seems 
to be homogeneous. The generalized canonical trans¬ 
formations eliminate the squeezing at each plane z and 
none signature of it is obtained. These insights lead us 
to consider the comoving frame as the physical one and 
therefore the proper one for quantization. Furthermore, 
interestingly, if we rewrite this Momentum in the original 
phase space, we obtain: 

Ma{q,P,z) = ^[{PaPa - qa PaY + PoaiQa/Paf]- (18) 

This is the space-analog of an Ermakov-Lewis invariant 
d]. So, it can be said that the Momentum (HH) is a spatial 
Ermakov-Lewis invariant in a comoving frame and that 
this invariant is the generator of spatial propagation in 
this kind of media. 

Now, following the principle of quantization of quan¬ 
tum mechanics (Qa,Pa) —t {Qa,Pa), where [Qa,P(T'] = 
ih5a,a', we obtain the Momentum operator: 

M, = \[Pl{s)+Pl,Ql{s)]. (19) 


In order to obtain the eigenvalues and eigenfunctions of 
this operator, we define the following annihilation and 
creation operators: 

^cr('S) = —7====[[JoaQ a + *Ar], (20) 

V ztipoa 

Atis) = - * 4 ]. ( 21 ) 

Therefore, equation (fT^ in terms of ([21]) can be rewritten 
as: 

j0l, = hPoa[AlA^ + l/2]. ( 22 ) 

The eigenstates of this Momentum are the same as those 
of the Number operator N^j = Aj.Ac, which fulfills the 
next eigenvalue equation: 

N^\N^) = N,\N,), (23) 

where and are eigenvalues and eigenstates for 
the Number operator, respectively. In terms of quantum- 
optical units | 20 j . we can redefine the quantum comov¬ 
ing scaled position (5 <t as the OFS operator whose 
eigenstates are those measured in homodyne experiments 
21 L l 22 j| and central to generalized quantum polarization 
that is £a = \/Poal^hQa- The eigenstates \Na) 
in the optical field-strength £ basis are given in terms of 
Hermite-Gauss functions as: 

21/4 

and their quantum propagation is given by a spatial 
Schrbdinger equation: 

jCia = flPoa [—7 2 ^a\ = f^Poa N ■ (25) 

‘*o£cr 

The solution of this equation yields the following propa¬ 
gated eigenstates \Na)'. 

01/4. i{N+l/2)6^ 

(26) 

with do- = Poa-Scr the total accumulated phase, in the 
same way as a quantum state propagating in a homo¬ 
geneous medium, but in this case in a comoving frame, 
which is the physical one for this kind of problems as 
shown above. Because of the similarity of this phase with 
the classical Gouy’s phase [1^, we call it spatial quantum 
Gouy’s phase. 

IV. PROPAGATION OF QUANTUM LIGHT IN 
LONGITUDINALLY INHOMOGENEOUS MEDIA 

Now, in order to gain insight into the behaviour of the 
longitudinally inhomogeneous waveguide, we study the 
spatial propagation of gaussian states, like the coherent 
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and squeezed quantum states of light, in the OFS rep¬ 
resentation. From the wavefunction (EH) and Mehler’s 
formula [ij, we easily obtain the propagator of the sys¬ 
tem for every mode [l4l| : 

K{S, £o; 0) = n (sin 

TT 

a 

sin&cr {S^+SQ^) — 2Ef7E0f7] (27) 

and the wavefunction of any pure state at the end of the 
medium can be worked out by means of: 

4'(£; z) = j K{e, So-, z, 0) 0) dSo, (28) 

where Sq and S stand for the A^-mode field-strength at 
the beginning and any point of the z-dependent medium, 
respectively. 



FIG. 1. Cosine-type effective refractive index N[z) = j5{z)/ko 
(solid line, left) and total accumulated classical phase 6 (dash 
line, right) versus propagation distance z, for parameters 
Nt = Pt/ko = 1.515, An = 0.5 and A = fco/50 nm~^ for 
a wavelength A = 653 nm. 


A single-mode gaussian state ja), where a = lale®*^, in 
the OFS space is given by [l^ : 


4'(£’o;0) = 


— to — a cost 


■AS, 


ASi 


-}e 


—iSo 


(29) 

where do = sint)) [|ap cos 2|q;|£’o], (f) = —ujt is the phase 
associated to the temporal evolution of the state at every 
plane z and ASq is the quantum noise or squeezing factor, 
with a value of one for a coherent state or different from 
one for squeezed states. Inserting equations EH and eh 
in (E51) . after a long but straightforward calculation, we 
obtain the wavefunction at any point z of the waveguide: 


^(S-,z) = ( 


^1/4 gie/2g-Z(5 


TT AS^ ' 

exp{-^[S-\a\ cos(<))-h6»)]^}, 


(30) 


where we have denoted 9,y = 9 and defined 5 = sin((f) + 
9) [jap cos{4> + 0) — 2|q;|£ 1] and: 

A£ = |^+cos=9Ar>|'y 

tand 

H = arctany j. 

Therefore, a coherent state keeps its quantum noise con¬ 
stant through propagation, AS'^ = ASq = 1, and gets 
a Gouy’s phase equal to the classical accumulated phase 
given by (IH, 0 = ^, as expected, unlike what we would 
obtain via the usual approach. On the other hand, an 
input squeezed state is both amplitude and phase modu¬ 
lated by the inhomogeneous medium, leading to different 
squeezing oscillations from those corresponding to an ho¬ 
mogeneous medium, and with a squeezing factor given by 
m- It is important to outline that the phase acquired 
by a squeezed state (l32l) is not predicted by the classical 
theory. 


(31) 

(32) 
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FIG. 2. Evolution of the quantum noise A£^ in a cosine-type 
longitudinally inhomogeneous medium for coherent (solid 
line) and squeezed (dash and dash-dot lines) states of light. 



Next, we show a specific example of the effects this 
kind of medium produce over a single-mode gaussian 
quantum state like the above introduced. We present a 
longitudinally inhomogeneous cosine-type medium with 
h{z) = cos(Az) and length L = 211 j A, as that depicted 
in Figure 1. The solution of equation m is given in 
terms of linearly independent Mathieu functions u{z) and 
v{z) [ 2 ^. Applying these solutions into equations (0 and 
EH, we can obtain numerically the value of 9 and, there¬ 
fore, the squeezing factor EH and the Gouy’s phase (15^ , 
which characterize the output quantum state after prop¬ 
agation in this medium. In Figures 2 and 3 we show 
the behaviour of the Gouy’s phase 0 and the quantum 
noise AS^, respectively, for three different input quan¬ 
tum states: the solid line stands for the propagation of a 
coherent state where ASq = 1 and the dash and dash-dot 
lines for squeezed states with ASq =3/2 and ASq =1/2, 
respectively. As can be checked in Figure 2, the planes 
where maximum or minimum squeezing of the quantum 
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noise are produced depend on the accumulated phase 9, 
being different from those corresponding to an homoge¬ 
neous medium, and the value of the input quantum noise 
A£q. Likewise, as was outlined above, the classical equa¬ 
tions do not predict the behaviour of states far from the 
classical, since squeezed and coherent states get differ¬ 
ent amounts of Guoy’s phase 0, as is sketched in Figure 
3. The classical phase 9 is recovered at the planes of 
maximum and minimum squeezing. Finally, we stress 
the interest this phase may have in quantum interferom¬ 
etry, since with an appropriate LI medium in one arm 
of an interferometer, quantum interference can be ad¬ 
justed controlling the strength of the quantum noise in a 
squeezed input state via equation (1321) . 



V. SUMMARY 


We have studied the propagation of quantum states 
of light in separable longitudinally inhomogeneous me¬ 
dia. We have proposed canonical transformations in a 
comoving frame and derived the generator of propaga¬ 
tion, the classical Momentum, realizing that it is a spa¬ 
tial Lewis-Ermakov invariant. We have quantized it and 
obtained the propagator in the OFS representation which 
is physically consistent. Furthermore, we have presented 
the propagation of a gaussian quantum state of light in 
a cosine-type refractive index medium and shown the 
net effect produced by this family of media: a quantum 
Gouy’s phase dependent on the input quantum state with 
application in quantum interferometry. 


FIG. 3. Evolution of the quantum Gouy’s phase 0 in a 
cosine-type longitudinally inhomogeneous medium for coher¬ 
ent (solid line) and squeezed (dash and dash-dot lines) states BIBLIOGRAPHY 

of light. 


[1] Yeon K. H., Kim H. J., Urn Gh. L, George T. F. and 
Pandey L. N. 1994 Phys. Rev. A 50 (2) 1035 

[2] Fernandez Guasti, M. and Moya-Cessa, H. 2003 J. Phys. 
A: Math. Gen. 36 2069 

[3] Pedrosa, LA. and Rosas, A. 2009 Phys. Rev. Lett. 103 
010402 

[4] Choi, J-R. 2004 Int. Journ. Theor. Phys. 43 (10) 2113 

[5] Lewis, H.R. 1967 Phys. Rev. Lett. 18 (13) 510 

[6] Krivoshlykov, S.G. and Sauter, E.G. 1992 Appl. Opt. 31 
(12) 2017 

[7] Moya-Cessa, H., Fernandez Guasti, M., Arrizon, V.M. 
and Chavez-Cerda, S. 2009 Opt. Lett. 34 (9) 1459 

[8] Luks A. and Perinova V. 2002 Progress in Optics 43 295 

[9] Linares J., Nistal M. C. and Barral D. 2008 New J. Phys. 
10 063023 

[10] Abram, 1. 1987 Phys. Rev. A 35 (11) 4661 

[11] Glauber, R.J. and Lewenstein, M. 1991 Phys. Rev. A 43 
(1) 467 

[12] Sodha M.S., Ghatak A.K. 1977 Inhomogeneous Optical 
Waveguides, Plenum Press NY 

[13] Marcuse, D. 1974 Theory of Dielectric Optical Waveg¬ 
uides, Academic Press, London 


[14] Linares J., Barral D. and Nistal M. C. 2012 J. Nonlin. 
Opt. Phys. Mat. 21 1250032 

[15] Pinney, E. 1950 Proc. Am. Mathem. Soc. 1 (5) 681 

[16] Mandel, L. and Wolf. E. 1995 Optical Coherence and 
Quantum Optics, Cambridge University Press, Cam¬ 
bridge 

[17] Chetouani, L., Guechi, L. and Hamman T.F. 1989 Phys. 
Rev. A 40 (3) 1157 

[18] Misner, C.W., Thorne, K.S. and Wheeler, J.A. 1973 
Gravitation, W.H. Freeman and Co., San Francisco 

[19] Takagi, S. 1990 Prog. Theor. Phys. 84 (6) 1019 

[20] Loudon, R. 1982 The Quantum Theory of Light, Claren¬ 
don Press, Oxford 

[21] Vogel, W. 1990 Phys. Rev. A 42 (9) 5754 

[22] Schleich, W.P. 2001 Quantum Optics in Phase Space, 
Wiley-VCH, Berlin 

[23] Linares, J., Barral D., Nistal M. C. and Moreno, V. 2011 
J. Mod. Opt. 58 711 

[24] Barral D., Linares J. and Nistal M. C. 2013 J. Mod. Opt. 
60 (12) 941 

[25] Boyd, R.W. 1980 J. Opt. Soe. Am. 70 (7) 877 

[26] Casperson, L.W. 1985 Appl. Opt. 24 (24) 4395 





